Si-compatible candidates for high-K dielectrics with the Pbnm perovskite structure 
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We analyze both experimentally (where possible) and theoretically from first-principles the di- 
electric tensor components and crystal structure of five classes of Pbnm perovskites. All of these 
materials are believed to be stable on silicon and are therefore promising candidates for high-if 
dielectrics. We also analyze the structure of these materials with various simple models, decompose 
the lattice contribution to the dielectric tensor into force constant matrix eigenmode contributions, 
explore a peculiar correlation between structural and dielectric anisotropies in these compounds 
and give phonon frequencies and infrared activities of those modes that are infrared-active. We find 
that CaZrC>3, SrZrC>3, LaHo03, and LaY03 are among the most promising candidates for high- if 
dielectrics among the compounds we considered. 

PACS numbers: 77.22.-d,77.55.df,85.50.-n 
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I. INTRODUCTION 

As a result of the ongoing down-scaling of complemen- 
tary metal-oxide-semiconductor (CMOS) integrated cir- 
cuits, the Si02 gate oxide of field effect transistors is 
getting thinner and thinner in every new generation of 
devices. 1 Therefore the leakage current due to quantum- 
mechanical tunneling through the dielectric interface is 
increasing. One way to reduce this current is to replace 
Si02 with a material that has a higher dielectric constant. 
Such a high- if dielectric layer with the same effective di- 
electric thickness (i.e., providing the same capacitance) 
could be physically thicker and thus reduce the gate leak- 
age. 

In order for this replacement material to be useful in 
practical applications on silicon, it also needs to be sta- 
ble in contact with silicon up to ~1000 °C, and among 
other things it must also have an appropriate band align- 
ment with silicon. 2 ~ 4 Currently, a hafnia-based dielec- 
tric is used as a replacement to Si0 2 in advanced CMOS 
transistors in production. 5 ~ 7 There are, however, draw- 
backs to this material too, e.g., the limited K that it 
provides and undesirable threshold voltage shifts arising 
from highly mobile oxygen vacancies. 8 This brings up the 
natural question: which other materials exist that would 
satisfy these requirements and would enable the scaling 
of MOSFETs to continue beyond today's hafnia-based di- 
electrics? 

The stability of single component oxides on silicon has 
been demonstrated both experimentally and from ther- 
modynamic analysis, 3 and a candidate list of multicom- 
ponent oxide materials has been compiled. 2 A promising 
group of these materials consists of perovskite oxides hav- 
ing a Pbnm (or closely related P2i/c) space group. These 
compounds are at the focus of the present work. Some 



of them have been studied in thin-film form, 9-27 but the 
full dielectric tensor of these materials has not yet been 
established, making the selection of materials best suited 
for high-if. applications difficult. Some of these mate- 
rials could also be of interest for microwave dielectric 
applications. 28,29 Thus we decided to study, both theo- 
retically and experimentally, the structural and dielectric 
properties of these compounds. The calculations are car- 
ried out using density-functional theory, and we compare 
the results with experimental data where we could obtain 
suitable samples for measurements. To our knowledge, 
previous theoretical calculations have been carried out 
in only a few cases. 30 ~ 32 

The paper is organized as follows. Explanations of 
both the experimental and theoretical methods used in 
this work are given in Sec. II. The main results on the 
structural and dielectric properties are given and dis- 
cussed in Sec. III. There we also discuss the correla- 
tions between the structural and dielectric properties of 
these perovskites, decompose the ionic contribution of 
the dielectric tensor into components arising from vari- 
ous force constant matrix eigenvectors, and discuss the 
effect of Ba antisite defects on the dielectric properties. 
We finish with a brief summary in Sec. IV. 

Supplementary material 33 contains the results of our 
calculations of the zone-center phonon frequencies, as 
well as the infrared activities for those modes that are 
infrared-active. 



A. Compounds under consideration 

In this work, we consider the following five groups of 
perovskites having the Pbnm space group. 

The first group are rare-earth scandates having for- 




(f) 




FIG. 1: Photographs of single crystal rare-earth scandates grown along the [110] direction by the Czochralski method, (a) 
PrScOa with a diameter of 12 mm, (b) NdSc03 with a diameter of 16 mm, (c) SmScOs with a diameter of 18 mm, (d) GdScC>3 
with a diameter of 32 mm, (e) TbScOs with a diameter of 18 mm, and (f) DySc03 with a diameter of 32 mm. 



mula ^4Sc0 3 where A is a rare-earth atom. In Sec. Ill C 1 
we report experimental measurements of the full dielec- 
tric tensors for PrSc0 3 , NdSc0 3 , SmSc0 3 , GdSc0 3 and 
DySc0 3 . Calculations were done on these and also 
on LaSc0 3 and TbSc0 3 . Note that HoSc0 3 , ErSc0 3 , 
TmSc0 3 , YbSc0 3 , LuSc0 3 , and YSc0 3 do not form sin- 
gle crystals with the perovskite structure from the melt at 
atmospheric pressure. Rather, they form solid solutions 
of A2O3 and Sc20 3 , i.e., (4,Sc) 2 3 , with the bixbyitc 
structure. 34 " 36 Nevertheless, LuSc0 3 27 and YbSc0 3 37 
have been formed in perovskite form as thin films via 
epitaxial stabilization, and others might be made in the 
same way. To analyze trends within this group of com- 
pounds, we also did the calculations of dielectric tensors 
on LuSc0 3 and YSc0 3 in the Pbnm perovskite structure; 
sec Sec. IIIC 1 for the details. 



The second group consists of rare-earth yttrates with 
formula AYO3. Only one such compound, LaY0 3 , is 
known to form a perovskite, 38 but to analyze trends the 
dielectric tensor of DyY0 3 in the perovskite structure 
was also calculated. 



In the third group we consider CaZr0 3 , SrZr0 3 and 
SrHf0 3 perovskites. Experimentally, we find that SrZr0 3 
and SrHf0 3 do not form single crystals, but instead are 
rather heavily twinned. 



The fourth group of compounds have the formula 
L^BB'Oe where the B atom is cither Mg or Ca and B' 
is either Zr or Hf. Little is known experimentally about 
these compounds, and single crystals of these compounds 
have not been made. 39 



The last group of compounds we considered have 
the formula AA'03, where both A and A' are rare- 
earth atoms. These include the 11 of such com- 
pounds that are known to form the perovskite struc- 
ture with space group Pbnm at atmospheric pres- 
sure: LaHo0 3 , LaEr0 3 , LaTm0 3 , LaYb0 3 , LaLu0 3 , 
CcTm0 3 , CeYb0 3 , CeLu0 3 , PrYb0 3 , PrLu0 3 , and 
NdLu0 3 . 38,40 ' 41 We calculated the dielectric and struc- 
tural properties of all of these compounds. The exper- 
imental determination, however, of the dielectric tensor 
in this group of compounds was done only for LaLu0 3 ; 
the results will be published elsewhere. 42 



II. PRELIMINARIES 
A. Structure of Pbnm perovskites 

The ideal cubic perovskite ^4i?0 3 consists of a network 
of corner-shared octahedra, each with an oxygen on its 
vertices and a B atom at its center, and A ions that are 
12-fold coordinated in the spaces between octahedra. It is 
well known that perovskites having sufficiently small A- 
site ions (i.e., a small Goldschmidt tolerance factor 43,44 ) 
often allow for a distorted perovskite structure that has 
a rotated framework of oxygen octahedra and displaced 
A-site ions. This lowers the space group symmetry from 
cubic Pm3m (0\) to orthorhombic Pbnm (D^), and 
the number of ABO3 formula units per primitive cell in- 
creases from Z = 1 to Z = 4, as shown in Fig. 2. The 
rotations of the octahedra in the Pbnm space group can 
be decomposed into two steps. The first step is the rota- 
tion around the [110] direction of the original cubic frame 
(the cubic frame is rotated by 45° around the z axis with 
respect to the Pbnm frame) by an angle 6> R , as in Fig. 3(a), 
and the second step is a rotation around [001] by 9m as 
in Fig. 3(b). The rotations must be done in that order 
to prevent distortions of the octahedra. The pattern of 
neighboring octahedral rotations is denoted by (a~a~c + ) 
in Glazer notation 45 (or see the directions of the arrows 
in Fig. 3). These rotations also allow for the displacement 
of ^4-site ions in the x-y plane without further lowering 
of the space-group symmetry. 



B. Computational methods 

The main computational method we are using is 
the density-functional theory as implemented in the 
Quantum-Espresso package. 46 The exchange-correlation 
functional was approximated using a generalized gradient 
approximation (GGA) of the Perdew-Burke-Ernzcrhof 
type 47 and ultrasoft pscudopotcntials were employed. 48 
The electronic wavefunctions were expanded in a basis of 
plane waves with kinetic energy up to 40 Ry, while the 
charge density was expanded up to 300 Ry. The Brillouin 
zone was sampled using a 4 x 4 x 4 Monkhorst-Pack grid. 49 

A new set of ultrasoft pseudopotentials 48 for the lan- 
thanoide series of rare earths, from La to Lu, were gener- 
ated for the present project. In all cases the /-shell filling 
was chosen as appropriate for the 3+ valence state: one / 
electron for Ce, two for Pr, etc. The / electrons were then 




FIG. 2: (Color online) 20-atom primitive cell of a Pbnm- 
distorted j4_B03 perovskite. A-site atoms are shown in black, 
B-site atoms in blue (at the centers of the octahedra) and 
oxygen atoms in red (at the vertices of the octahedra). Or- 
thorhombic unit cell vectors (a, b, and c) are also indicated. 
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FIG. 3: (Color online) (a) Projection on the a-c plane of the 
structure with Ok > and 9m — 0. (b) Projection on the a-b 
plane for (9r, = and 9m > 0. Color coding of atoms and 
positions of axes labels are the same as in Fig. 2. In (a), a 
darker shading is used to indicate the two octahedra that are 
further away along the y-coordinate; in (b), the two bottom 
octahedra are exactly below the two top octahedra. 



considered to be in the core (and un-spin-polarizcd) for 
the proposes of generating the pseudopotentials. Thus, 
the / electrons are not explicitly included in the solid- 
state calculations. Such an approximation can be jus- 
tified whenever the strong on-site Coulomb interactions 
of electrons in the / shell drive the occupied / states 
well below, and the unoccupied states well above, the 
energy range of interest for spd bonding in the crystal. 
Of course, this will not be a good approximation for some 
heavy-fermion or mixed-valent systems, and in any case 
our approach is obviously unable to describe phenomena 
involving magnetic ordering of / electrons at low tem- 
perature. Nevertheless, we believe that this approach is 
quite reasonable for the present purposes. 

The artificial nature of the scattering in the / channel 



did, however, pose some problems in the pscudopotential 
construction. In particular, we found that the lattice con- 
stant of a perovskite containing the rare-earth atom in 
question could differ for two pseudopotentials having dif- 
ferent scattering properties in the / channel; this causes 
problems since the usual approach of matching to the 
all-electron / scattering is not appropriate in the present 
case. To ameliorate this problem, the /-channel parame- 
ters of these pseudopotentials were optimized so that re- 
sulting pseudopotentials would give the "correct" cell vol- 
umes for simple rare-earth compounds. Since the GGA 
typically overestimates crystal volumes by about 1-2%, 50 
the optimization was actually done in order to produce 
a corresponding overestimate in a consistent fashion. 

For this procedure, our compounds of choice were the 
rare-earth nitrides with the simple rock-salt structure. 
The experimental 51 and calculated volumes of these ni- 
trides are indicated in Fig. 4. Note that the volumes 
of CeN, PrN, and GdN show an anomalous behavior 
that is presumably due to strong correlation effects as- 
sociated with the proximity to a mixed-valent regime, 52 
and therefore they will not be correctly treated by our 
GGA calculation. To avoid this problem we first carried 
out a smoothened fit of the experimental volumes versus 
atomic number over the lanthanoide nitride series, but 
with CeN, PrN, and GdN omitted from the fit, as shown 
by the solid line in Fig. 4. We then used these fitted val- 
ues to set the target volumes for the optimization of the 
pseudopotentials. 
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FIG. 4: Unit cell volumes of all rock-salt rare-earth nitrides, 
in A . Empty squares are the results obtained using our op- 
timized ultrasoft pseudopotentials for rare-earth atoms; solid 
circles are experimental 51 results. The solid line is a fit to the 
experimental values excluding CeN, PrN, and GdN. 

We used density-functional perturbation theory 53 to 
calculate the dielectric response. Both purely electronic 
e cl and ionic e Ion contributions were calculated. 54 The 
electronic part is defined as 
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where P a is the polarization induced by the electric field 



£g while all ions are held fixed (u = 0). The remaining 
component of the dielectric response is by definition the 
ionic contribution e lon . 

This ionic part can be calculated from the force- 
constant matrix §i a ,jp and the Born effective charge ma- 
trix Zi. a p. The force-constant matrix is defined as 
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where E is the total energy of the system and Ui a is the 
displacement of the «-th atom along the direction a. We 
will denote the n-th normalized eigenvector of this matrix 
as g? a and its eigenvalue as /i„ . The Born effective charge 
matrix is defined as 
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where P a is the polarization induced in a crystal by the 
displacement of the i-th atom in the direction j3. V is the 
volume of the unit cell and e is the electron charge. Fi- 
nally, the ionic part of the dielectric tensor can be written 

as 
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where the charge Q™ of the n-th eigenmode is defined 



through the effective charge matrix as Q™ 



i(i 



Zi,ap£™p- 



C. Experimental methods 

1. Crystal growth 

PrSc0 3 , NdSc0 3 , SmSc0 3 , GdSc0 3 , TbSc0 3 , and 
DySc0 3 single crystals were grown using an auto- 
mated Czochralski technique with RF-induction heat- 
ing. 55 ' 56 Pre-dried powders of Sc20 3 , Pr 6 On, Nd20 3 , 
Sni2 3 , Gd2 3 , Tb407, and Dy20 3 , were mixed in 
the stoichiometric ratio, pressed, and sintered at about 
1400 °C for 15 h. Due to the high melting temperature 
of PrSc0 3 , NdSc0 3 , SmSc0 3 , GdSc0 3 , TbSc0 3 , and 
DySc0 3 (~2100 °C), a crucible (cylindrical with 40 mm 
or 60 mm diameter and 40 mm or 60 mm height, depend- 
ing on the crystal diameter) and an active afterheater 
made of iridium were used. Flowing nitrogen or argon 
was used as the growth atmosphere. Those of these rare- 
earth scandates for which the radiative heat transport 
via the crystal is hindered by absorption suffer from a 
serious problem of bulk crystal growth, namely that they 
tend to exhibit a spiral growth which distinctly decreases 
the yield. 55 Due to the lack of seed crystals, the initial 
growth experiments were performed with an iridium seed 
rod. Because these materials tend to grow as large single- 
crystalline grains, suitable seeds could be selected at a 
very early stage. 



All rare-earth scandatc crystals were grown along the 
[110] orientation. The pulling rate was 0.8-2 mm-h -1 
and the rotation rate was 8-15 min -1 (depending on the 
crystal diameter). The crystals were 35-50 mm in length 
and 18 or 32 mm in diameter. The PrSc0 3 crystals 
have a green color, the NdSc0 3 crystals are dark purple, 
the SmSc0 3 and DySc0 3 crystals are light yellow, the 
GdSc0 3 crystals arc colorless, and the TbSc0 3 crystals 
are nearly colorless (see Fig. 1). 



In addition to the Czochralski growth of PrSc0 3 , 
NdSc0 3 , SmSc0 3 , GdSc0 3 , TbSc0 3 , and DySc0 3 , the 
single-crystal growth of HoSc0 3 , YSc0 3 , and solid so- 
lutions of rare-earth scandates that would be expected 
to have smaller lattice constants than DySc0 3 (e.g., 
Dyo.5Luo.5Sc0 3 ), was also attempted using the floating- 
zone technique. In all cases, pre-dried powders of the 
rare-earth oxides (R2O3) and Sc20 3 were ground, pressed 
into a rod, and sintered at 1400 °C for 6 hours. The 
floating-zone machine was a Gero (model SPO), and 
heating was performed with two bulbs at the focal point 
of ellipsoidal mirrors. Unfortunately, the 1 kW rated 
output power of the quartz-halogen lamps used in this 
floating-zone system was close to the melting temper- 
ature of the compounds, limiting the size of the crys- 
tals that could be grown to < 2 mm diameter. SrZr0 3 , 
SrHf0 3 , and LaSc0 3 single crystals grown by the floating 
zone technique were also studied. 57 



The available i?20 3 -Sc20 3 binary phase diagrams 
show a distinct transition between the Dy20 3 -Sc20 3 and 
the Ho20 3 -Sc20 3 systems (Dy and Ho are neighbor- 
ing elements in the periodic table). According to pub- 
lished phase diagrams, the rare-earth elements of La, 
Nd, Sm, Gd, and Dy all form i?Sc0 3 compounds that 
melt congruently. In contrast the oxides of Ho, Er, Tm, 
Yb, and Lu form complete solid solutions with Sc20 3 
at their melting points. 34 ~ 36 (The same happens with 
Y as well. 34-36 ) Although HoSc0 3 has been synthesized 
with solid-state techniques at temperatures (well below 
the melting point) at which the perovskitc polymorph of 
HoSc0 3 is stable, our results confirm that the perovskite 
polymorph of HoSc0 3 is not stable at its melting point 
in agreement with existing phase diagrams. Analysis of 
our Y-Sc-0 and Ho-Sc-0 single crystals revealed that 
they were solid solution mixtures of Ho20 3 and Sc20 3 , 
and Y 2 3 and Sc 2 3 , not the desired perovskite poly- 
morphs of HoSc0 3 and YSc0 3 . The only compound 
that did not behave according to its phase diagram was 
LaSc0 3 . Attempts to grow LaSc0 3 resulted in crystals 
that were a mixture of three different phases including 
LaSc0 3 , Sc20 3 , and a third, unidentified phase. 57 X- 
ray analysis of the post-annealed polycrystalline feed rod 
showed single-phase LaSc0 3 , consistent with the existing 
phase diagram. The reason for the introduction of sec- 
ond phases during melting is not clear and requires more 
study. 



2. Electrical characterization 

Because PrScOg, NdSc0 3 , SmScOg, GdScOg, and 
DyScOa arc orthorhombic at room temperature, their 
dielectric tensor contains three independent coefficients 
that can be measured along the three principle crys- 
tal axes. 58 From the grown PrScC>3, NdSc0 3 , SmSc0 3 , 
GdScC>3, and DyScC>3 single crystals, slices were cut 
in different orientations for electrical characterization. 
For GdSc03, samples were cut from the as-grown sin- 
gle crystals along the [100], [010], and [001] directions, 
and the dielectric tensor coefficients were measured us- 
ing a parallel-plate capacitor configuration. For PrScC>3, 
NdSc03, SmScOs, and DyScOs, several additional ori- 
entations were prepared, so that a least-squares fit tech- 
nique could be applied using standard methods 58 to cal- 
culate the dielectric tensor, using their known lattice 
parameters. 56 ' 59 

Unfortunately, the SrZrC>3 and SrHfOs crystals were 
heavily twinned. The formation of twins in the three 
pscudocubic orientations is facilitated in these two or- 
thorhombic systems because the deviation from the cu- 
bic symmetry is very small. As a result, only an average 
dielectric constant could be measured. 

To obtain capacitors, gold or platinum electrodes were 
evaporated onto both sides of the approximately 0.5 mm 
thick slabs with areas ranging from 15 to 100 mm 2 , and 
capacitance measurements were made at room temper- 
ature with an HP4284A using a 16034E test fixture at 
10 kHz, 100 kHz, and 1 MHz. No edge capacitance cor- 
rections were used because of the large ratio of the elec- 
trode area to the thickness of the samples. The room 
temperature dielectric loss of the samples was gener- 
ally very small (< 0.1% for most samples) and the fre- 
quency dispersion was negligible in the measured range. 
The temperature dependent dielectric measurements on 
GdScO-3 and DySc0 3 were made using a HP 4284a LCR 
meter with a dipstick cryostat. 



III. RESULTS AND DISCUSSION 



2 • 10~ 3 . The volume of the unit cell, on the other hand, 
is consistently overestimated by 1-2%, as is usually ex- 
pected from the GGA exchange-correlation functional, 
and as we would expect from our construction of the 
rare-earth pseudopotentials. 

All structures show an angle 6*r that is about y/2 times 
larger than 9m- Therefore, consecutive rotations by #r 
and Oyi can be considered approximately as a single ro- 
tation around a {111} axis in the cubic frame. That is, 
the actual (a~a~c + ) pattern of rotations is very nearly 
(a~a~a + ) in the Glazer notation. 45 See Sec. IIIF for a 
more detailed discussion. 



1. Rare-earth scandates 

The rare-earth scandates ^4Sc03 show a decrease in 
volume by ~9% while going along the series from A=L& 
to yl=Dy (the calculated primitive unit cell volume is 
271.40 A 3 for LaSc0 3 and 249.81 A 3 for DySc0 3 ). On 
the other hand, the Sc-0 distance remains nearly con- 
stant along the series (2.12 A for LaSc03 and 2.11 A 
for DyScOs), which means that the change in volume 
is almost entirely due to the larger octahedral rotation 
angles for DyScOs as compared to LaScOs. Our calcu- 
lations also show that the same trend continues all the 
way to LuSc03. 



2. Rare- earth yttrates 

The rare-earth yttrates have a very similar behavior as 
the rare-earth scandates. The main quantitative struc- 
tural difference between the two comes from the fact that 
yttrium is a larger ion than scandium. This leads to a 
larger volume for the yttrates, and also a larger rotation 
angle due to a smaller tolerance factor. 



A. Structural properties 

We focus first on the structural properties of these sys- 
tems. The structure of the Pbnm perovkites is described 
by three orthorhombic lattice constants plus two A-site 
and five oxygen Wyckoff parameters. Figure 5 shows 
graphically the most important structural parameters of 
these systems, while Table I gives detailed information 
on all the structural parameters. Rotational angles in 
Table I and in Fig. 5 were calculated by fitting the struc- 
tural parameters to a model in which the octahedra are 
perfectly rigid (see Sec. IIIB for the details of this model). 

Overall we find good agreement with experimental val- 
ues for the structural parameters. The Wyckoff coordi- 
nates in particular are in excellent agreement with ex- 
periments, with the average error being on the order of 



3. CaZrOs, SrZrOs, and SrHfO-j, 

SrZrOs and SrHfOs have quite similar structural prop- 
erties. The main difference can be traced to the fact that 
Hf is a smaller ion than Zr. Therefore, the calculated 
average Hf-0 distance is 2.07 A, while the average Zr-0 
distance is 2.11 A. Furthermore, their octahedral rotation 
angles are about 1.7 times smaller than in the rare-earth 
scandates. 

In CaZrOs the average Zr-0 distance is 2.10 A, which 
is very close to the corresponding distance in SrZr03 
and SrHfOs. Thus, the main reason why CaZr03 has 
a smaller volume than SrZrOs is because of the larger 
rotation angles in CaZrOs. 



4- La-zBTi Og compounds 



9m (see Fig. 3) leads to Wyckoff parameters given by 



1 



We consider L^BB'Oq compounds with B—Mg or Ca 
and _B'=Zr or Hf. These compounds are expected to ex- 
hibit rock-salt ordering of the B-site ions as a result of 
the difference in charge and ionic radius between the B 
and B' ions. 60 This ordering reduces the symmetry from 
the orthorhombic Pbnra to the monoclinic P2\/c (C| ;i ) 
space group. 

The structural properties for these systems are re- 
ported in Fig. 5 and in Table II. The rotational angles are 
obtained by a fit to the rigid-octahedra model in which 
we have allowed for different sizes of B- and ^-centered 
octahedra. (See the end of Sec. IIIB for details.) 

The unit cell volume is larger by about 5 A 3 per prim- 
itive cell for the compounds containing Zr than for those 
containing Hf. On the other hand, compounds with Ca 
are larger by about 28 A 3 than those containing Mg. 
Similarly, the rotation angles are larger in compounds 
containing Ca than in those with Mg. The discrepancy 
between octahedral sizes is largest for La2CaHfOg (12% 
linear increase) and smallest for I^MgZrOg (0.4% linear 
increase). 



X2 = 



tan#R, 



5. Rare-earth rare-earth perovskites 



We now briefly analyze the structural properties of 
Pbnm perovskites of type AA'O?, where both A and A" 
are rare-earth atoms. All eleven compounds we consid- 
ered are known experimentally to form the pcrovskitc 
structure in the Pbnm space group. 38,40,41 

Among these 11 compounds, the largest unit-cell vol- 
ume of 311.58 A 3 is found in LaHo03, and the small- 
est of 292.32 A 3 is in NdLu03, Oxygen oxtahcdral rota- 
tion angles are quite large in all of these compounds and 
show very little variation from one compound to another. 
The trends of the rotation angles are as expected from a 
tolerance-factor analysis: perovskites with smaller ^4-site 
ions but the same 5-site ions have larger oxygen octa- 
hedral rotation angles, and the opposite is true for the 
5-sitc ions. 



B. Comparison with model of perfectly rigid 
octahedra 
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Here we have denoted the Wyckoff coordinates of the 
oxygen atoms at the 4c Wyckoff point with xi and y%, 
while those of the remaining oxygen atoms at the 8d point 
are denoted with X3, 2/3, and 23. The Wyckoff coordinates 
of the ^4-site ion at the 4c point arc denoted by X\ and 
2/1, but these are left unspecified in our rigid-octahedra 
model. It also leads to orthorhombic lattice constants 
given by 




COS#rCOS#m, 
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(12) 

where Vq is the volume the structure would have if the 
octahedra were rotated rigidly back to #r = 9m = 0. 

The Wyckoff parameters and unit-cell ratios from our 
calculations can be well fitted by Eqs. 5-12 (see Table I 
and Fig. 5 for the values of the fitted angles). By far 
the largest discrepancy is found for Wyckoff parameter 
g/2- For a typical system (e.g., LaScOa) the discrepancy 
between calculated and fitted 2/2 values is about 0.016, or 
50% with respect to the difference from the cubic case. 
For the remaining oxygen Wyckoff coefficients, the dis- 
crepancy averages about 0.003, or ~5%. 

The rotation angles for the La2-RB'06 systems were ob- 
tained by fitting their structural parameters to a slightly 
more complicated model of rigid octahedra than the one 
given in Eqs. 5-12. In this model, we first change the 
relative sizes of B- and -B'-centered octahedra. The ratio 
of their linear sizes is denoted by d/df. We then proceed 
with the rotation by an angle #r. around the [110] axis 
in the cubic frame. Finally, we perform a rotation of the 
-B-centcred octahedra around [001] by an angle 9m, and 
of the -B'-centered octahedra by an angle 9' M around the 
same axis. The resulting fitted values of these parameters 
are given in Table II. 



C. Dielectric properties 



In Pbnm perovskites, a rigid rotation of the oxygen 
octahedra by (9r followed by another rigid rotation by 



In this section we discuss the dielectric properties of the 
materials included in our study. Pbnm perovskites are 



TABLE I: Structural parameters of the Pbnm perovskites we considered. Our calculated values are indicated with the letter T, 
experimental values with E and other theoretical data with 0. Fitted octahedra rotation angles #r and #m are given in degrees, 

see Sec. IIIB for the details. A-site ions occupy the 4c site and in cubic configuration xi — and y\ — 1/2. _B-site ions are on 

the 4a site. One type of oxygen sites are at 4c and in the cubic case X2 — and j/2 = 0. The remaining oxygens are at 8d and 

X3 — 1/4, yz = 1/4 and zs — in the cubic case. 

Unit cell parameters Wyckoff coordinates Model 

a b c V b/a c/a xi yi X2 2/2 X3 J/3 23 Or 8m 

(A) (A) (A) (A a ) (°) (°) 

LaSc0 3 T 5.7030 5.8414 8.1469 271.40 1.0243 1.4285 0.0123 0.5467 -0.0977 -0.0323 0.2059 0.2943 0.0523 14.9 9.9 

E a 5.6803 5.7907 8.0945 266.25 1.0194 1.4250 0.0100 0.5428 -0.0968 -0.0277 0.2073 0.2958 0.0521 14.6 9.6 

PrScOs T 5.6372 5.8367 8.0908 266.21 1.0354 1.4352 0.0148 0.5537 -0.1049 -0.0377 0.2014 0.2979 0.0557 16.4 11.0 

E j 5.608 5.780 8.025 260.1 1.0307 1.4310 0.0121 0.5507 -0.1052 -0.0395 0.1977 0.3008 0.0555 16.2 11.2 

NdScOs T 5.6077 5.8317 8.0667 263.80 1.0399 1.4385 0.0159 0.5562 -0.1083 -0.0401 0.1997 0.2992 0.0574 17.0 11.5 

E b 5.577 5.777 8.005 257.9 1.0359 1.4354 0.0133 0.5532 -0.1088 -0.0418 0.1953 0.3020 0.0571 16.8 11.8 

SmSc0 3 T 5.5483 5.8067 8.0196 258.37 1.0466 1.4454 0.0176 0.5596 -0.1159 -0.0455 0.1964 0.3019 0.0610 18.2 12.5 

E b 5.531 5.758 7.975 254.0 1.0410 1.4419 0.0149 0.5566 -0.1163 -0.0468 0.1935 0.3037 0.0609 17.9 12.5 

GdSc0 3 T 5.4987 5.7794 7.9861 253.79 1.0510 1.4524 0.0191 0.5617 -0.1222 -0.0502 0.1941 0.3036 0.0642 19.0 13.3 

E b 5.481 5.745 7.929 249.7 1.0482 1.4466 0.0163 0.5599 -0.1209 -0.0501 0.1912 0.3052 0.0628 18.7 13.2 

TbScOs T 5.4764 5.7646 7.9735 251.72 1.0526 1.4560 0.0198 0.5624 -0.1251 -0.0524 0.1932 0.3043 0.0656 19.4 13.7 

E k 5.4543 5.7233 7.9147 247.07 1.0493 1.4511 0.0167 0.5603 -0.1239 -0.0545 0.1900 0.3054 0.0643 19.1 13.7 

DyScOg T 5.4560 5.7501 7.9629 249.81 1.0539 1.4595 0.0203 0.5630 -0.1276 -0.0545 0.1923 0.3050 0.0669 19.7 14.1 

O f 5.449 5.739 7.929 248.0 1.0532 1.4551 0.019 0.562 -0.130 -0.057 0.190 0.307 0.068 19.9 14.1 

E b 5.443 5.717 7.901 245.9 1.0503 1.4516 0.0174 0.5616 -0.1262 -0.0561 0.1886 0.3063 0.0659 19.4 13.9 

LaY0 3 T 5.9035 6.1225 8.5810 310.16 1.0371 1.4535 0.0173 0.5506 -0.1284 -0.0581 0.1948 0.3035 0.0689 18.9 13.6 

E 1 5.890 6.086 8.511 305.1 1.0333 1.4450 

CaZrOs T 5.5974 5.7875 8.0416 260.51 1.0340 1.4367 0.0133 0.5506 -0.1078 -0.0413 0.1976 0.2999 0.0572 16.6 11.6 

E m 5.5831 5.7590 8.0070 257.45 1.0315 1.4341 0.0122 0.5495 -0.1044 -0.0401 0.1976 0.3000 0.0554 16.2 11.4 

SrZrOs T 5.8068 5.8602 8.2323 280.14 1.0092 1.4177 0.0070 0.5311 -0.0759 -0.0195 0.2140 0.2857 0.0401 11.6 7.7 

O s 5.652 5.664 7.995 255.9 1.0021 1.4145 0.007 0.534 -0.107 -0.036 0.199 0.301 0.056 14.7 10.0 

E c 5.7963 5.8171 8.2048 276.65 1.0036 1.4155 0.0040 0.5242 -0.0687 -0.0133 0.2154 0.2837 0.0363 10.4 7.2 

SrHfOs T 5.7552 5.7754 8.1365 270.45 1.0035 1.4138 0.0052 0.5230 -0.0660 -0.0128 0.2209 0.2792 0.0346 10.0 6.2 

O h 5.6887 5.7016 8.0455 260.95 1.0023 1.4143 0.006 0.528 -0.073 -0.016 0.2166 0.2834 0.0385 10.8 7.0 

E d 5.7516 5.7646 8.1344 269.70 1.0023 1.4143 0.0040 0.5160 -0.0630 -0.0140 0.2189 0.2789 0.0335 9.6 6.4 

LaHoOa T 5.9135 6.1367 8.5859 311.58 1.0377 1.4519 0.0170 0.5508 -0.1293 -0.0587 0.1936 0.3044 0.0692 19.0 13.6 

E> 5.884 6.094 8.508 305.1 1.0357 1.4460 

LaEr0 3 T 5.8971 6.1174 8.5509 308.48 1.0374 1.4500 0.0169 0.5509 -0.1272 -0.0567 0.1945 0.3036 0.0681 18.7 13.3 

E 1 5.870 6.073 8.465 301.8 1.0346 1.4421 

LaTmOs T 5.8829 6.1002 8.5190 305.73 1.0369 1.4481 0.0167 0.5509 -0.1252 -0.0549 0.1951 0.3030 0.0670 18.5 13.1 

E 1 5.859 6.047 8.453 299.5 1.0321 1.4427 

LaYbOs T 5.8692 6.0833 8.4890 303.09 1.0365 1.4464 0.0164 0.5509 -0.1233 -0.0532 0.1957 0.3024 0.0659 18.3 12.8 

E 1 5.843 6.033 8.432 297.2 1.0325 1.4431 

LaLuOs T 5.8579 6.0695 8.4646 300.95 1.0361 1.4450 0.0162 0.5508 -0.1217 -0.0518 0.1962 0.3020 0.0651 18.2 12.6 

E c 5.8259 6.0218 8.3804 294.00 1.0336 1.4385 0.0138 0.5507 -0.121 -0.056 0.193 0.307 0.063 17.9 12.7 

CeTm0 3 T 5.8520 6.0870 8.4984 302.72 1.0401 1.4522 0.0179 0.5528 -0.1283 -0.0575 0.1942 0.3036 0.0686 19.0 13.5 

E 1 5.828 6.035 8.405 295.6 1.0355 1.4422 

CeYbOa T 5.8381 6.0707 8.4671 300.09 1.0398 1.4503 0.0176 0.5529 -0.1264 -0.0558 0.1948 0.3031 0.0676 18.8 13.3 

E 1 5.806 6.009 8.373 292.1 1.0350 1.4421 

CeLu0 3 T 5.8270 6.0578 8.4420 297.99 1.0396 1.4488 0.0174 0.5529 -0.1249 -0.0544 0.1952 0.3027 0.0667 18.7 13.1 

E 1 5.793 5.997 8.344 289.9 1.0352 1.4404 

PrYbOa T 5.8085 6.0544 8.4481 297.10 1.0423 1.4544 0.0185 0.5542 -0.1296 -0.0585 0.1936 0.3039 0.0692 19.3 13.7 

E 1 5.776 5.995 8.368 289.8 1.0379 1.4488 

PrLuOa T 5.7974 6.0424 8.4217 295.01 1.0423 1.4527 0.0183 0.5543 -0.1281 -0.0571 0.1936 0.3036 0.0683 19.1 13.5 

E 1 5.768 5.991 8.340 288.2 1.0387 1.4459 

NdLuOa T 5.7699 6.0270 8.4062 292.32 1.0446 1.4569 0.0193 0.5555 -0.1310 -0.0596 0.1932 0.3042 0.0699 19.5 14.0 

E 1 5.737 5.974 8.311 284.8 1.0413 1.4487 

a Reference [61]. b Reference [62]. c Reference [63]. d Reference [64]. ° Reference [40]. f Reference [30]. s Reference 

[31]. h Reference [32]. 'Reference [38]. j Reference [59]. k Reference [65]. 'Reference [66]. m Reference [67]. 



TABLE II: Structural parameters of the L^BB'Oq perovskites we considered, space group P2i/c. Angles are given in degrees, 
unit cell vectors in A and unit cell volumes in A . Monoclinic angle between a and c lattice vectors is denoted as a ac - j4-site 
ions occupy the 4e site with coordinate (3/4, 1/2, 1/4). B-site ions (either Mg or Ca) occupy the 2a site and B' (Zt or Hf) 
occupy the 2d site. There are three non-equivalent positions for oxygens and they all occupy 4e sites. Coordinates of oxygen 
atoms in the cubic case are (1/2, 0, 1/4), (1/4, 1/4, 0), and (3/4, 1/4, 0). Fitting parameters Or, Om, 0' m and d/d' are also given; 
see Sec. IIIB for the details. 
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FIG. 5: Structural information for all systems we considered. Bottom pane shows volume in A 3 per formula unit (f. u.) of 
j4-B03. Calculated values are shown as red circles, and experimental values as blue plus symbols, if available. Top pane shows 
oxygen octahedra rotation angles in degrees. Theoretical values are shown with red circles and squares and experimental values 
with blue plus and cross symbols. Or angles are shown with circles and plus symbols while <9m angles with squares and cross 
symbols. For 'LslzBB'Os systems the average of 8m and 0' M is given. For numerical values see Table I and II. 



orthorhombic and thus have diagonal dielectric tensors, 
with e xx y^ e yy 7^ e zz in general. In addition to reporting 
these components, we also focus on analyzing the results 
in terms of the three linear combinations 
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representing the average dielectric tensor, a measure of 
the x-y anisotropy, and a measure of z anisotropy, re- 
spectively. This choice of parameters was made to sim- 
plify the analysis of trends of dielectric properties of 
these compounds. The theoretical - and where available, 
experimental - results for the dielectric-tensor compo- 
nents are reported in Figure 7 and in Table III. The 
theoretical values are further decomposed in Table III 
into purely electronic or frozen-ion contributions e cl and 
lattice- mediated contributions e lon . We find that the elec- 
tronic contribution is roughly five times smaller than the 
ionic one, is nearly isotropic, and does not show a dra- 
matic variation from one perovskite to another. Thus, 
it is clear that the lattice-mediated ionic contributions 
play by far the dominant role in the observed dielectric 
tensors and their anisotropics. 

Our calculations of the zone-center phonon frequencies 
as well as the infrared activities for those modes that 
are infrared-active are given in the supplementary mate- 
rial. 33 

We now consider each of our chosen classes of Pbnm 
perovskitcs in turn, orienting the presentation from the 
point of view of the theoretical calculations, but mention- 
ing the comparison with experiment where appropriate. 



1. Rare-earth scandates 

All rare-earth scandates ^4ScC>3 have rather similar val- 
ues for their isotropically-averaged dielectric constants, 
falling between about e = 26 and e = 28. The xx compo- 
nent for all these systems is larger than the yy component 
by about Aey = 4. On the other hand, the zz compo- 
nent changes significantly from LaSc03 to DySc03. In 
LaSc03 the average of the xx and yy components is al- 
most as large as the zz component (Aej_=-1), while in 
DyScOs, the zz component is larger by about Ae^=9 
than the average of xx and yy components. 

These results are in good agreement with experiment, 
especially for e and Ae||. On the other hand, Ae^ is 
consistently larger in experiments by about 3-5, but the 
trend of increasing Ae±_ is present in both theory and 
experiment. 

As was mentioned earlier, rare-earth atoms heavier 
than Dy (i.e., Ho-Lu) and Y itself do not form single- 
crystal scandates. Nevertheless, at least some (YbSc03 37 
and LuSc03 27 ) can form Pbnm perovskites in thin- film 
form. In order to establish the trends of the dielectric 
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FIG. 6: Measured temperature dependence at 10 kHz of 
the dielectric tensor components (top panel) and dielectric 
loss (bottom panel) of GdSc03 (blue) and DyScOs (red). 
Downward-pointing triangles denote e xx and tanS xx ; upward- 
pointing triangles denote e yy and t&nS yy ; circles denote e 22 
and t&n8 zz . 



properties for these materials, we calculated the dielec- 
tric tensors of L11SCO3 and YSCO3. The dielectric ten- 
sor of LuSc03 shows the continuation of the trend from 
LaScOs to DyScOs. Both xx and yy components are 
slightly smaller than for DyScOs, their numerical values 
being 23.5 and 21.4 respectively. On the other hand, the 
zz component (44.8) is larger than for DyScOs (32.6) 
and for LaScOs (27.4). YSCO3 has dielectric tensor com- 
ponents of 26.9, 23.0, and 37.7 for its xx, yy, and zz 
components, respectively. 

The experimentally measured dependence of the di- 
electric tensor (and loss) components on temperature 
is shown in Fig. 6 for two compounds, GdSc03 and 
DyScOs. In both cases we find that the dielectric tensor 
properties do not change significantly with temperature 
over the examined range (4.2-470 K). The e zz compo- 
nent slightly decreases with temperature while the e xx 
and e vv components show the opposite behavior. Wc 
expect that similar trends will be observed in all other 
rare-earth scandates. 

As our calculated dielectric tensor is at K and our 
measured dielectric tensor (Table III and IV) is at room 
temperature, the absence of a significant temperature de- 
pendence of the dielectric tensor is important to our abil- 
ity to make a meaningful comparison between the calcu- 
lated and measured coefficients. 

Anomalies were observed in the dielectric loss along 
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the z axis (tan<5 22 ) for both DySc03 and GdSc03 below 
50 K. The origin of these is unknown, but their pres- 
ence suggests the possibility of a low-temperature phase 
transition. 



2. 



earth yttrates 



those with Mg, and a bit larger for those with Zr than for 
those with Hf. Therefore, the dielectric response in this 
class of materials is largest for La2CaZrC>6, with e = 28.5, 
and smallest for La2MgHfOg, with e = 23.6. All com- 
puted dielectric tensor components for these systems are 
given in Table IV. 



We now consider the rare-earth yttrates, i.e., ^4Y03 
where A is one of the rare-earth atoms. These are simi- 
lar to the rare-earth scandates, but with yttrium on the 
B site instead of scandium. Only one such compound, 
LaYC>3, is known to form a pcrovskitc, 38 but others 
might form in thin films via epitaxial stabilization. We 
find that LaY03 has a larger zz component than does 
LaScC>3. In LaYC>3 the zz component of the dielectric 
tensor is 38.0, while in LaScC>3 it is 27.4. On the other 
hand, the xx and yy components are almost unchanged 
with respect to LaSc03. In LaY03 the xx component 
is 30.6 and the yy component is 25.7. We also find that 
heavier rare-earth atoms on the j4-site tend to destabi- 
lize this Pbnm structure even further. For example, we 
find that DyY03 in the Pbnm structure has an unstable 
mode at z70 cm -1 that is IR-active along the z direction. 



5. Rare-earth rare- earth perovskites 

The 11 rare-earth-rare-earth perovskites we consid- 
ered show a bigger variation in the isotropically-averagcd 
dielectric constant e than do the rare-earth scandates 
(LaSc03 - DyScOs). The largest average dielectric con- 
stant among them is 32.9 in LaHo03. The largest com- 
ponent of a dielectric tensor is also found in LaHo03, 
whose e 22 is 41.7. 

The measure Aen of x-y anisotropy shows little varia- 
tion among the components in this series. The anisotropy 
is of the same sign as in the rare-earth scandates. 

Finally, the z anisotropy Aej_ once more shows 
a larger variation than in the rare-earth scandates. 
This anisotropy is largest for LaHo03 and smallest for 
LaLu03. 



3. CaZrOz, SrZr0 3 , and SrHfOs 

According to our calculations, SrZr03 and SrHfOs 
show large and rather isotropic dielectric tensors. The 
average dielectric tensor e is 40.9 and 32.8 in SrZrOs 
and SrHfOs, respectively. Their x-y anisotropics have 
an opposite sign as compared to all of the other com- 
pounds we analyzed. Unfortunately, because of twinning 
(see section II C 2), we could only measure an average 
dielectric constant for these systems, and therefore we 
could not directly compare our full calculated dielectric 
tensors with experiment. Still, if we make a comparison 
between theory and experiment for the average dielectric 
tensor e, the agreement is reasonable. 

Our calculations suggest that the CaZr03 compound, 
on the other hand, has a very high value of the z- 
anisotropy of 28.1. Its x-y anisotropy of 1.1, on the 
other hand, is quite small. The average dielectric ten- 
sor e = 43.6 is the highest among the all compounds we 
considered, mostly because of the very large e 22 compo- 
nent of the dielectric tensor. Very similar results were 
also obtained in other theoretical studies. 68 ' 69 



4- La-^BE! Og compounds 

The La2-RB'06 systems show a small, non-zero off- 
diagonal e xz component, —0.4 for La2MgZr03 and 4 for 
La2CaZr03. e xz is allowed because the space group is re- 
duced from orthorhombic {Pbnm) to monoclinic (P2i/c) 
for these compounds. Their isotropically-averagcd dielec- 
tric tensors are larger for systems containing Ca than for 



D. Decomposition of the ionic contribution to the 
dielectric tensor 

As already mentioned, the ionic contribution to the di- 
electric tensor dominates in all of the systems we consid- 
ered. The expression for the ionic contribution given in 
Eq. (4) provides a decomposition into contributions com- 
ing from eigenmodes of the force-constant matrix. The 
Pbnm symmetry in perovskites, which is also approxi- 
mately satisfied in La2BB'Oe compounds, allows a given 
cigenmodc to contribute only to a single component (e„, 
e yy , or e zz ) of the dielectric tensor. This decomposition 
is given in Fig. 8 for all three components. 

In rare-earth scandates, all three directions are evi- 
dently very different. The e xx component is dominated 
by a low-lying mode whose contribution is almost con- 
stant along the series (it contributes to e xx by 9.6 for 
LaSc03 and 11.1 for DyScOs). The e yy component, on 
the other hand, has sizable contributions coming from 
several modes. Finally, the e 22 component comes mostly 
from a single low-lying mode. Unlike for the e xx com- 
ponent, the contribution from the mode responsible for 
the e zz component changes dramatically across the series, 
varying from 6.9 for LaScOs to 16.3 for DyScOs. This 
explains the large value of the z anisotropy in DyScOs 
as compared to LaSc03 that is visible in Fig. 7. 

A behavior similar to that of the rare-earth scandates 
is also observed in LaY03 and in the rare-earth rare- 
earth perovskites. The SrZr03 and SrHf03 compounds 
show a quite similar behavior to each other. The e xx 
component has contributions coming from many modes, 
the e yv component is dominated by a single low-lying 
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TABLE III: Dielectric parameters of the perovskites with space group Pbnm that we considered. Our calculated values are 
denoted by the letter T, our experimental values with E, and other theoretical data with 0. First all three non-zero dielectric 
constant tensor components are given. Next, the average dielectric constant tensor, the x-y anisotropy, and z anisotropy are 
given; see Eqs. 13-15. Finally, the electronic and ionic contributions are given separately. 
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LaSc0 3 


T 


30.4 


26.4 


27.4 


28.1 


4.0 


-1.0 


5.0 


5.0 


4.8 


25.3 


21.5 


22.6 


PrSc0 3 


T 


28.6 


24.2 


26.2 


26.3 


4.4 


-0.1 


5.0 


5.0 


4.8 


23.5 


19.2 


21.5 




E» 


25.4 


27.3 


29.6 


27.4 


-1.9 


3.3 














NdScOs 


T 


27.8 


23.4 


26.1 


25.7 


4.4 


0.5 


5.0 


4.9 


4.7 


22.8 


18.5 


21.3 




E b 


25.5 


21.5 


26.9 


24.6 


4.0 


3.4 














SmSc0 3 


T 


27.2 


22.9 


27.3 


25.8 


4.3 


2.3 


4.9 


4.9 


4.7 


22.2 


18.0 


22.7 




E'' 


23.1 


19.9 


29.0 


24.0 


3.2 


7.5 














GdSc0 3 


T 


26.4 


22.5 


29.3 


26.1 


3.9 


4.9 


4.9 


4.8 


4.0 


21.6 


17.7 


24.7 




E'' 


22.8 


19.2 


29.5 


23.8 


3.0 


8.5 














TbSc0 3 


T 


26.1 


22.4 


30.7 


26.4 


3.7 


0.5 


4.8 


4.8 


4.0 


21.2 


17.6 


26.1 


DyScOs 


T 


25.7 


22.3 


32.6 


26.9 


3.5 


8.0 


4.8 


4.8 


4.5 


20.9 


17.5 


28.0 




O" 


24.1 


21.2 


27.7 


24.3 


2.9 


5.1 


4.9 


4.9 


4.7 


19.2 


16.3 


23.0 




E'' 


21.9 


18.9 


33.8 


24.9 


3.0 


13.4 














LaY0 3 


T 


30.6 


25.7 


38.0 


31.4 


4.9 


9.9 


4.7 


4.0 


4.3 


25.9 


21.1 


33.7 


CaZrOs 


T 

E' 


34.8 


33.7 


62.4 


43.6 
30 


1.1 


28.1 


4.0 


4.7 


4.0 


30.2 


29.1 


57.8 


SrZrOa 


T 


38.0 


41.5 


43.3 


40.9 


-3.4 


3.5 


4.0 


4.0 


4.0 


33.4 


36.9 


38.7 




O' 


19.9 


21.5 


23.0 


21.5 


-1.6 


2.3 


5.1 


4.9 


4.8 


14.8 


16.6 


18.2 




E b 








32 


















SrHf0 3 


T 


30.0 


35.1 


33.2 


32.8 


-5.1 


0.7 


4.3 


4.3 


4.3 


25.7 


30.9 


29.0 




O'' 


33.1 


46.8 


40.8 


40.2 


-13.7 


0.9 


4.4 


4.3 


4.3 


28.7 


42.4 


36.5 




E b 








26.2 


















LaHo0 3 


T 


31.0 


26.0 


41.7 


32.9 


5.0 


13.3 


1.7 


4.0 


4.3 


26.3 


21.4 


37.4 


LaEr0 3 


T 


29.9 


25.0 


36.4 


30.4 


4.9 


9.0 


4.7 


4.0 


4.3 


25.2 


20.4 


32.1 


LaTm0 3 


T 


29.0 


24.2 


33.3 


28.8 


4.8 


0.7 


4.0 


4.0 


4.3 


24.4 


19.7 


29.0 


LaYbOs 


T 


28.3 


23.6 


30.9 


27.6 


4.7 


5.0 


4.0 


4.5 


4.3 


23.6 


19.1 


26.6 


LaLuOs" 


T 


27.6 


23.1 


29.3 


26.7 


4.5 


3.9 


4.0 


4.5 


4.3 


23.0 


18.6 


25.0 


CeTm0 3 


T 


27.9 


23.5 


34.9 


28.8 


4.4 


9.2 


4.0 


4.0 


4.3 


23.3 


18.9 


30.6 


CeYb0 3 


T 


27.2 


22.8 


32.1 


27.4 


4.4 


7.0 


4.0 


4.0 


4.3 


22.6 


18.3 


27.8 


CeLu0 3 


T 


26.6 


22.3 


30.2 


26.4 


4.3 


5.7 


4.0 


4.5 


4.3 


22.0 


17.8 


25.9 


PrYb0 3 


T 


26.6 


22.7 


34.7 


28.0 


4.0 


10.1 


4.0 


4.5 


4.3 


22.0 


18.2 


30.4 


PrLuOs 


T 


26.1 


22.0 


32.2 


26.8 


4.0 


8.1 


4.0 


4.5 


4.3 


21.5 


17.5 


27.9 


NdLuOs 


T 


25.3 


21.7 


34.5 


27.2 


3.0 


11.0 


4.0 


4.5 


4.3 


20.8 


17.2 


30.2 


"Reference [30]. 
6 This work. 
c Reference [31]. 
Reference [32]. 
^Experimental data 
^Reference [67]. 


in reference [42]. 























TABLE IV: Dielectric parameters of the La2-B-B'Oe perovskites we considered. First all four dielectric constant tensor compo- 
nents are given, followed by the average value of the dielectric constant tensor. Finally, the electronic and ionic contributions 
to dielectric constant tensor are given, x axis is chosen along lattice vector a, y along b, and z close to c (with the small 
component along a due to monoclinic cell). 
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4.6 
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-0.06 
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18.9 


19.8 


-0.3 


La 2 MgHf0 6 


24.9 


22.9 


22.9 


-0.3 


23.6 


4.6 


4.5 


4.4 


-0.05 


20.3 


18.4 


18.5 


-0.3 


La2CaZrOe 


29.9 


24.9 


30.6 


3.8 


28.5 


4.7 


4.6 


4.4 


-0.06 


25.2 


20.3 


26.2 


3.8 


La 2 CaHf0 6 


27.4 


23.0 


26.6 


2.0 


25.7 
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4.3 


-0.04 
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18.5 


22.3 
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FIG. 7: Dielectric information for all systems we considered. Bottom pane shows average dielectric tensor e, middle pane shows 
x-y anisotropy of the dielectric tensor, Aey , and z anisotropy of the dielectric tensor Aex is given in top pane. Calculated 
values are shown with red circles and experimental values with blue cross symbols, if available. See Table III for numerical 
values. 



mode, and e zz is dominated by two low-lying modes. On 
the other hand, e zz component in CaZr03 shows very 
large contribution coming from a single low-lying mode. 
Finally, we note that the La2BB'Oe compounds contain- 
ing Ca have stronger contributions to e xx and e zz from 
low-lying modes than do those containing Mg. 



E. Compounds with R3c symmetry 

At room temperature the ground state of BiFeOs is 
ferroelectric with polar space group R3c, the pattern of 
octahedral rotations being (a~a~a~) in the Glazer no- 
tation. At higher temperature, however, BiFeC-3 under- 
goes a phase transition in which the fcrroclcctricity and 
the (a~a~a~) pattern of octahedral rotations disappear 
simultaneously. 70,71 This observations led us to hypoth- 
esize that rotations of octahedra around the pseudocu- 
bic [111] axis, as in the (a~a~a~) pattern, tend to be 
energetically compatible with the presence of a ferroelec- 
tric distortion along the same axis. This would tend to 
suggest that perovskitcs that adopt the centrosymmctric 
R3c group, which also exhibits the (a~a~a~) pattern of 
oxygen octahedra, might be close to a ferroelectric insta- 



bility leading to the lower-symmetry R3c space group, 
and thus that such compounds might have an especially 
large component of the dielectric tensor along the pseu- 
docubic [111] axis. 

To test this hypothesis, we have carried out a scries 
of calculations on SrZrC>3 and GdScC>3 in which struc- 
tural relaxation was allowed while maintaining the R3c 
symmetry. In both compounds we find some IR-active 
phonon modes that either have very low or imaginary 
frequency, indicating a near or actual instability. In the 
case of SrZrC>3 we find a mode that is active along the 
[111] pscudocubic direction and has an extremely small 
frequency of only 6 cm -1 , while for GdSc03 we find that 
the corresponding mode is unstable with an imaginary 
frequency of il42 cm -1 . These calculations show that 
imposing the R3c structure on SrZrOs and GdScC>3 make 
them nearly or actually ferroelectric, thus confirming our 
hypothesis. 

Incidentally, the observation that SrZr03 is more likely 
than GdSc03 to be stabilized in the R3c structure is 
consistent with the fact that perovskite structures that 
prefer smaller rotation angles are more likely to form R3c 
than Pbnm structures, as discussed by Woodward 44 . We 
find that the rotational angles for SrZrOs in the Pbnm 
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space group are t9pt = 11.6° and 9m = 7.7°, while in 
GdScOs they are substantially larger, 6*r, = 19.0° and 
9m = 13.3°. More directly, we also find that the ground- 
state energy of SrZr0 3 having the i?3c structure is only 
higher by 33.6 meV per formula unit than in the Pbnm 
structure. On the other hand, in GdScC>3 the R3c is 
higher in energy by a much larger increment of 386 meV. 
Finally, we note that LaA10 3 , NdA10 3 , and BaTb0 3 
may also be of interest, as these all have the i?3c space- 
group symmetry and should also be chemically stable on 
silicon. 



F. Correlation between structural and dielectric 
properties 

The heuristic observation about BiFe0 3 mentioned in 
the previous section (Sec. HIE) led us to make a more 
detailed analysis of the correlation between structural 
and dielectric anisotropics in all five groups of Pbnm 
perovskites. As can be seen from Fig. 3, the presence 
of the octahedral rotations breaks the symmetry among 
the three Cartesian directions in the Pbnm perovskites. 
One would therefore naively expect that the anisotropy 
in the dielectric tensor component should also be corre- 
lated with the size of these rotation angles, but this is 
not what we observe. For example, LaSc0 3 and DySc0 3 
both have rather substantial octahedral rotation angles 
(0 R is 14.9° in LaSc0 3 and 19.7° in DySc0 3 ), but they 
have very different values of the dielectric z anisotropy 



(Ae^ is —1.0 in LaSc0 3 and 8.6 in DySc0 3 ). An even 
more extreme behavior can be seen in the case of rare- 
earth perovskites between, e.g., LaHo0 3 and LaLu0 3 . 

Thus, wc find no simple correlation between the dielec- 
tric tensor anisotropies and the values of the octahedral 
rotation angles. Instead, we find a correlation between 
the dielectric tensor anisotropies and the mismatch of the 
two rotation angles 9r and 9m-, as we explain next. 

While the Pbnm symmetry does not impose any re- 
lationship between the two octahedral rotation angles 
0r and 9m, we find in practice that all the compounds 
we studied obey the heuristic relationship 9r ~ v2$M- 
This means that the oxygen octahedra are rotated about 
the three Cartesian axes by almost the same rotation 
angle, or cquivalently, that the rotation axis is nearly 
(111). In the Glazer language, these Pbnm perovskites 



having (a 
to an (a 



a c + 

n /7~i~ 



) rotations can be said to be very close 
) pattern. We can measure the mis- 



match between the actual (a~a~a + ) and the hypothet- 
ical (a~a~c + ) rotation pattern by the quantity 9m — 
0r/v2, and it is this quantity that we find to be strongly 
correlated with the dielectric anisotropy Aej_ . 

This is shown in Fig. 9, where Aej_ is plotted ver- 
sus 9m — 9^/^/2 for all of the compounds considered 
in this work. It is apparent that the Ill-III-valent per- 
ovskites have a different behavior than the II-IV-valent 
ones. Nevertheless, we conclude that in both cases there 
is a strong correlation between the mismatch angle and 
the dielectric tensor anisotropy. The sign of the correla- 
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tion is such that a deviation from the (a~a~a + ) pattern 
having an increased rotation angle around the z axis gives 
a larger dielectric tensor component along the z axis, and 
thus a larger z anisotropy Aej_. 
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FIG. 9: Correlation between the dielectric tensor z-anisotropy 
(Aex) and the mismatch in the oxygen octahedra rotation an- 
gle. Only some perovskites are labeled. Perovskites where 
both A and B ions are III-valent are indicated with blue 
square symbols and those where A ions are II-valent and B 
are IV-valent are indicated with red circles. 



G. Antisite substitutions 

Experimentally the compositions of the perovskites 
that we have been describing up to now by their nominal 
compositions, e.g., LaLu03, are in fact slightly different 
from the compositions of the single crystals on which the 
dielectric tensors were measured. This is because our 
crystals are grown at the congrucntly melting composi- 
tions, e.g., Lao.94Lu1.06O3, which differ from the nomi- 
nal compositions described up to now. The congruently 
melting compositions of all relevant Pbnm perovskites 
studied have been found to be poor in the ^4-site cation 
and rich in the -B-site cation composition. 59 ' 72,73 

For this reason, we decided to carry out a theoretical 
analysis of the effects of B atoms substituting at the A 
site on the structural and dielectric properties of the ma- 
terial. Detailed calculations were done only for the case 
of LaLuC>3, but we expect that similar trends will be ob- 
served in the remaining rare-earth rare-earth perovskites 
as well as in the rare-earth scandates and yttrates. Of 
course, other kinds of compositional disorder might also 
be present, but such possibilities are not analyzed here. 



1. Analysis of antisite defects in LaLuOs 

We studied LuL a antisite defects in LaLuC>3 using a 
superccll approach. Specifically, in order to model a sit- 
uation in which one of every 16 La atoms is substituted 



by Lu, which is about a 6% substitution, we constructed 
an 80- atom supercell containing a single antisite defect. 
The supercell is enlarged with respect to the primitive 
20-atom Pbnm primitive cell by doubling along both the 
orthorhombic a and b lattice vectors. The resulting stoi- 
chiometry is 

(Lao.9375Luo.0625)Lu03 Or Lao.9375Lui.o62503. 



The presence of the LuL a antisite in this particular 
80-atom supercell reduces the crystal symmetry from or- 
thorhombic Pbnm to monoclinic Pm. After full relax- 
ation of the crystal structure in this space group, we find 
that the a, b, and c lattice vectors are reduced by 0.3%, 
0.2% ,and 0.1%, respectively, while the monoclinic angle 
between a and b lattice vectors of 90.03° deviates only 
very slightly from 90°. 

The influence of the Lul & substitution on the dielectric 
properties is more complex. Evaluated in the same co- 
ordinate frame as in the Pbnm unit cell, the e xx and e yy 
dielectric tensor components remain almost unchanged, 
and the new e xy component allowed by the monoclinic 
symmetry is quite small, only 2.1. On the other hand, the 
e zz component is drastically altered by the presence of Lu 
atom on the La site. In fact, we find that the 80-atom su- 
percell is actually just barely unstable in the Pm space 
group, as indicated by the presence of a phonon mode 
with a very small imaginary frequency of il6 cm -1 . 74 The 
contribution of this phonon mode to the e zz component 
(evaluated in the unstable Pm structure) is therefore neg- 
ative, specifically, —33.6. Since this phonon frequency is 
so close to zero, we expect that it would get renormalized 
to positive frequency at room temperature. For this rea- 
son, we did not follow the structural relaxation of our 80- 
atom supercell along the direction of the unstable mode, 
and a realistic estimate of the dielectric response of the 
system is difficult. Nevertheless, we conclude that LuLa 
substitutions in LaLuOs have the potential to increase 
the e zz dielectric tensor component substantially. 



2. Discussion 



As can be clearly seen in Fig. 7, our calculated z 
anisotropy (Aej_) is consistently larger than the mea- 
sured one for all the rare-earth scandates for which we 
have experimental measurements. In view of the calcula- 
tions reported for LaLuOs above, we tentatively attribute 
this discrepancy to the generic tendency of B atoms to 
substitute on the A site in these compounds. This ob- 
servation is consistent with the fact that smaller B ions 
that substitute for larger A ions will reside in a relatively 
larger cage, providing room to rattle and thereby con- 
tribute to an enhanced dielectric response. 
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IV. SUMMARY 

The main focus of this work has been the application of 
both computational and experimental methods to study 
the structural and dielectric properties of various Pbnm 
perovskites that have potentially large dielectric tensor 
components and are chemically stable on silicon up to 
^1000 °C. 2 Such compounds might be good candidates 
for future use as high- if dielectrics in microelectronics ap- 
plications, e.g., as a possible replacement of hafnia-based 
high- K dielectrics currently used in the CMOS transis- 
tors in integrated circuits. 

Of the compounds we have considered, CaZr03, 
SrZr03, LaHo03, and LaY03 appear to be especially 
promising. CaZrOs has the largest calculated average 
dielectric tensor (e = 43.6) among the compounds we 
considered, and SrZr03 is a close second with e = 40.9. 
The dielectric tensor in CaZr03 is very anisotropic, with 
its e zz component almost twice as large as e xx or e yy , 
while on the other hand SrZr03 has an almost isotropic 
dielectric tensor. Unfortunately, the full dielectric ten- 
sors of these compounds have not yet been measured due 
to lack of single crystals. 

Of the rare-earth rare-earth Pbnm perovskites, only 
LaLu03 has had its dielectric tensor measured to date. 
Our results on this compound will be presented in de- 
tail elsewhere. 42 The theoretical calculations, however, 
indicate that other compounds in this series should have 
even larger dielectric tensor components, with LaHo03, 
having e = 32.9, being the most promising among these. 
LaY03 is expected to behave very similar to LaHo03 
since Y and Ho have almost the same ionic radii, so it 
may be promising as well (e = 31.4). Thus, this series of 
compounds clearly deserves additional scrutiny. 

Of course, there are good reasons for preferring amor- 
phous over single-crystalline materials for such high-if 
applications. Certainly the ability of amorphous Si02 



to conform to the substrate and to eliminate electri- 
cal traps played a central role in its dominance as the 
gate dielectric of choice for 40 years for silicon-based 
metal-oxide-semiconductor ficld-cffcct transistors. The 
present hafnia-based high- if dielectrics are amorphous 
or nanocrystalline. 75 For this reason, any eventual appli- 
cation of these materials for high- if applications would 
presumably require the adoption of one of two strate- 
gics. The first is the possibility of growing crystalline 
epitaxial oxides directly on silicon, which clearly would 
require a very high level of control of interface chem- 
istry and morphology before it could become a practi- 
cal solution. The second is the possibility that some of 
the compounds investigated here could be synthesized in 
amorphous or nanocrystalline form. We have not inves- 
tigated these issues here, nor have we tried to calculate 
what (possibly very substantial) changes in the dielectric 
properties would occur in the amorphous counterparts, 
as this would take us far beyond the scope of the present 
study. Nevertheless, these are important questions for 
future investigations. 
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